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, Abstract. We report on the 20 ksec observation of Vela 
' X-1 performed by BeppoSAX on 1996 July 14 during its 
[ Science Verification Phase. We observed the source in two 
' intensity states, characterized by a change in luminosity of 
I a factor ~ 2, and a change in absorption of a factor ~ 10. 
' The single Narrow Field Instrument pulse-averaged spec- 
, tra are well fit by a power law with significantly different 

■ indices. This is in agreement with the observed changes of 
, slope in the wide-band spectrum: a first change of slope 

at ^ 10 keV, and a second one at ~ 35 keV. To mimic 
this behaviour we used a double power law modified by 
' an exponential cutoff — the so-called NPEX model — 
I to fit the whole 2-100 keV continuum. This functional is 
' able to adequately describe the data, expecially the low in- 
I tensity state. We found an absorption-like feature at ^ 57 

■ keV, very well visible in the ratio performed with the Crab 
] spectrum. We interpreted this feature as a cyclotron res- 
' onance, corresponding to a neutron star surface magnetic 

strength of 4.9 x 10^^ Gauss. The BeppoSAX data do not 
require the presence of a cyclotron resonance at ^ 27 keV 
as found in earlier works. 

Key words: binaries: eclipsing — Stars:individual 
(HD77581) — Stars:neutron — pulsars:individual (Vela 
X-1) — X-rays:stars 



1. Introduction 

Among the class of X-ray binary pulsars. Vela X-1 is 
the one that has been more carefully monitored since its 
discovery as a pulsator in 1976 (McClintock et al. 1976). 
It is the prototype of the subclass of pulsars that 
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accrete matter directly from the stellar wind, com- 
ing from the massive early-type star IID77581 (spec- 
tral type BO lb (Hutchings 1974)). During its orbit 
around the companion, the X-ray emission from the pul- 
sar suffers substantial photoelectric absorption, increas- 
ing until the source is eclipsed by its companion ev- 
ery 8.96 days — the orbital period — for about 1.7 
days (Avni 1976). The study of the orbital phase depen- 
dence of the low energy absorption due to the circum- 
stellar matter, together with spectroscopic studies of the 
optical companion, has shown the presence of accretion 
(Blondin, Stevens & Kallman 1991) and photoionization 
(Kaper, Hammerschlag-Hensberge & Zuiderwijk 1994) 
wakes. 

The pulse-averaged spectrum of Vela X-1 also 
shows an Iron emission line centered at ^ 6.4 keV 
(White, Swank & Holt 1983). Recently, cyclotron reso- 
nance features (CRFs) at ~ 27 and 55 keV been reported 
by Kretschmar et al. (1996) and Makishima & Mihara 
(1992). 

2. Observation 

The BeppoSAX satellite is a program of the Italian 
Space Agency (ASI), with participation of the Nether- 
lands Agency for Aerospace Programs (NIVR), devoted 
to X-ray astronomical observations in the broad 0.1- 
300 keV energy band (Boella et al. 1997a). The pay- 
load includes four Narrow Field Instruments (NFIs) and 
two Wide Field Cameras (Jager et al. 1997). The NFIs 
consist of Concentrators Spectrometers (C/S) with 3 
units (MECS) operating in the 1-10 keV energy band 
(Boella et al. 1997b) and 1 unit (LECS) operating in 0.1- 
10 keV (Parmar et al. 1997), a High Pressure Gas Scin- 
tillation Proportional Counter (HPGSPC) operating in 
the 3-120 keV energy band (Manzo et al. 1997) and a 
Phoswich Detection System (PDS) with four scintillation 
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detection units operating in the 15-300 keV energy band 
(Frontcra ct al. 1997). 

During the Science Verification Phase a series of well 
known X-ray sources have been observed in order to 
check the capabilities and performances of the instruments 
aboard BeppoSAX . Vela X-1 is one of these sources and 
it has been observed by three of the four NFIs (LEGS 
was not operative during this pointing) on 1996 July 14 
from 06:01 to 20:54 UT. The net exposure time for MEGS, 
HPGSPG and PDS was 21.6, 8.1 and 9.5 ksec, respectively. 
These differences arc due to rocking of collimators (sec be- 
low) and different filtering criteria during passages in the 
South Atlantic Geomagnetic Anomaly (SAGA) and before 
and after Earth occultations. 

In Fig. 1 we show the average light curve of Vela X-1 
as observed by the three operative NFIs instruments. Data 
were tclcmctrcd in direct modes, which provide complete 
information on time, energy and, if available, position for 
each photon. The two collimated instruments (HPGSPG 
and PDS) were operated in rocking mode with a 96 sec 
stay time for HPGSPG, and 50 sec stay time for PDS, in 
order to monitor the background along the orbit — a com- 
plete on-sourcc observation, apart from SAGA passages 
and Earth occultations, is possible only for PDS because 
of its two (mechanical) collimators. The single-collimator 
HPGSPG will always show stay time gaps in its data. 

The first part of the observation (about 40% of the to- 
tal net time, and marked LOW in the figure) corresponds 
to one of the common intensity dips showed by Vela X -1, 
explained as due to passages of the neutron star through 
clumpy circumstellar material (Sato et al. 1986). 

3. Spectral Analysis 

The pulse-averaged spectrum of Vela X-1 has been fitted 
separately in the energy bands of the three operative Bep- 
poSAX NFIs instruments, and the results are summarized 
in Table 1. We used XSPEC V9.0 (Arnaud 1996) as spec- 
tral fitting program, and we refer to its User Manual^ for 
the details of the functional used for the fits. We treated 
separately the two intensity states of the source. Note that 
the high xi for PDS spectra is mainly due to systematic 
uncertainty, below 5% in the fit residuals. 

From Table 1 it is evident that any attempt of fit- 
ting the wide-band Vela X-1 average spectrum with a 
single power law, eventually modified by energy cutoff 
and/or narrow-band absorption features, will be unsuc- 
cessful, because of the completely different spectral slopes 
in the three NFIs energy ranges. Indeed, the usual — 
for an X-ray pulsar — power law plus high-energy cutoff 
(White, Swank & Holt 1983) gives a very poor fit to the 
2-100 keV spectra, with xl ^ ^ and 9 for the two inten- 
sity states (Orlandini et al. 1997a). A fit with a different 
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^ The XSPEC User Manual is available on-line at 
http://legacy.gsfc.nasa.gov/docs/xanadu/xspec/u_manual.html 
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Time (sec since 14 Jul 1996 00:00 UT) 

Fig. 1. Vela X-1 background subtracted light curve for all the 
three operative BeppoSAX NFIs instruments. Gaps are due 
to South Atlantic Anomaly passages and Earth occultations 
(the average BeppoSAX duty cycle is about 50%). The upper 
scale represents the orbital phase based on the ephemeris given 
by Deeter et al. (1987). The first panel shows the light curve, 
rebinned at 283 sec, from one of the three MEGS aboard Bep- 
poSAX. The second panel shows the 384 sec HPGSPG hght 
curve, while in the third panel the 500 sec PDS data are dis- 
played. In both last cases binning times have been chosen as 
an integer multiple of the rocking stay time. 



continuum, namely a power law modified by the so called 
Fermi-Dirac cutoff (Tanaka 1986), gave an unacceptable 
fit, too. 

Therefore we used for the continuum a spectral law 
that allowed the flattening of the spectrum observed in 
the ^ 10-30 keV range. In particular we found that two 
power laws modified at high energy by the same exponen- 
tial cutoff were able to describe the 2-100 keV spectra, 
expecially the LOW state. This spectral functional has 
been chosen as a compromise between a reasonable good 
fit and a simple law. It is also the same continuum the 
so-called NPEX — used by Mihara in fitting the contin- 
uum of X-ray pulsars (Mihara 1995). 

By including a Gaussian line in emission at ^ 6.4 
keV and an absorption edge at ~ 7.7 keV we ob- 
tained a reasonable good fit, with x^= 1-3 (331 dof) 
for LOW state, and xl= 2.8 (332 dof) for HIGH 
state. From the form of the residuals at high en- 
ergy we were led to add GRFs in the spectrum 
description. Other authors (Makishima & Mihara 1992, 
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Table 1. BeppoSAX NFIs spectral fits to Vela X -1 pulse-averaged spectrum. MECS data are fit by a single power law, plus 
photoelectric absorption, an emission Gaussian line and an absorption edge. The HPGSPC spectrum is fit by a power law with 
exponential cutofT plus an emission Gaussian line (the inclusion of a CRF did not improve the fit). The PDS data have been fit 
with a power law modified by a high-energy cutoff of the form exp[(£c — E)/Ef], plus a Lorenzian CRF. 
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Note — All quoted errors represent 90% confidence level for a single parameter. 

" Total fiux in each instrument energy band, in units of photons cm~^ scc^'. 

'' Total X-ray luminosity in each instrument energy band, in units of 10^*" ergs sec^', assuming a distance of 1.9 kpc 
(Sadakane et al. 1985). 



Kretschmar et al. 1996) have considered the fundamen- 
tal cyclotron resonance at ^ 27 kcV, but our high en- 
ergy data from both HPGSPC and PDS do not show it. 
Therefore we added to our model a CRF at ~ 57. Both 
a Lorenzian (Mihara 1995) or a Gaussian in absorption 
(Soong et al. 1990) model gave a quite good fit, expe- 
cially for the LOW state. The fit results are summarized 
in Table 2 and the fits for the two intensity states are 
shown in Fig. 2. The inter-calibration between the NFIs 
is quite good, expecially at the light of this early calibra- 
tion status^. With MECS2 as reference, normalization is 
0.98 for HPGSPC, and 0.79 for PDS. 

Our results on the Iron line energy at 6.4 keV are in 
agreement with those found by other experiments. A more 

detailed discussion on the properties of the Iron line, with 
emphasis on its pulse-phase variation, will be discussed in 
a separate paper. 



Updated information on the calibration status are available 
at http://www.sdc.asi.it/software/ 



4. Discussion 

4-1. The X-ray continuum 

This is the first time that we have three (quite different) 
X-ray detectors that simultaneously observed Vela X-1 in 
a so wide energy range. From the single NFIs spectra we 
can extrapolate that the wide-band spectrum will have a 
complex form: up to ~ 10 keV it is very well described by a 
simple power law (MECS data). Then the spectrum begins 
to flatten: first at ~ 10 keV (cutoff in HPGSPC data) and 
then more pronouncedly after ~ 35 keV (cutoff in PDS 
data). These variations are quite in agreement with the 
change observed in the Vela X-1 pulse profile, from a five- 
peaked to a two-peaked shape just in the 10-30 keV range 
(see e.g., Orlandini 1993). 

The NPEX functional that we used in the wide-band 
spectral fitting is surely a crude approximation of this be- 
haviour, but it describes the data very well for the LOW 
state, and reasonably well for the HIGH state. It is impor- 
tant to stress how the choice of the functional describing 
the continuum has dramatic effects on the line character- 
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Fig. 2. BeppoSAX Vela X-1 pulse-averaged spectra in the two intensity states: LOW on the left and HIGH on the right. The 
continuum has been modelled by two power laws plus the same exponential cutoff. We added an emission line, an absorption 
edge and a Lorenzian cyclotron resonance at ~ 57 keV. Note the dramatic change in the low energy absorption. 



Table 2. Fit results on the wide-band BeppoSAX Vela X-1 
pulso-avcragcd spectra in the two intensity states. Also an ab- 
sorption Gaussian at 57±1 keV, with width 22±1 keV FWHM, 
and Equivalent Width 30 ±3 keV fits well the CRF. All quoted 
errors represent 90% confidence level for a single parameter. 
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~ 20 40 60 80 100 

Energy (keV) 

Fig. 3. Upper panel: Ratio between the PDS Vela X-1 and 
Crab spectra. Note the change of slope at ~ 57 keV. Middle 
panel: the Vela X 1/Crab ratio multiplied by E-"'^, the func- 
tional form of the Crab spectrum. Deviations from the con- 
tinuum (dotted line) occur for E 35 keV. Lower panel: ratio 
between the (Vela X-l/Crab ratio x E-"'^) and the functional 
form of the Vela X-1 continuum as derived from the wide-band 
fit. A Gaussian fit to the CRF (dotted line) is also shown. 



istics. Indeed, a fit to our data with a power law modified 
by a high energy cutoff of the type described by White, 
Swank & Holt (1983) was not able to well describe both 
the ~ 27 and ~ 57 keV CRFs (Orlandini et al. 1997a). 



4-2. Cyclotron resonance features 

The Vela X-1 pulse-averaged spectrum shows all the fea- 
tures present in a typical X-ray binary pulsar spectrum, 
namely an Iron emission line with its absorption edge, 
and signatures of cyclotron resonance. While the Iron line 
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feature is clearly determined in the spectrum, there are 

ambiguities in the position of the fundamental CRF. In- 
deed, a fit to the wide-band BeppoSAX spectrum with an 
NPEX and two CRFs, with the fundamental constrained 
in the 10-40 keV range, produced a cyclotron energy at 
~ 24 keV and a xt comparable with those shown in Ta- 
ble 2, although the CRF depth was extremely low. An 
F-test was not successful in discriminating between the 
two hypothesis, with a probability of chance improvement 
of 98% for LOW State and 42% for HIGH State. 

But there arc some arguments against the low energy 
CRF: first, the HPGSPC would surely be successful in its 
detection if it had the same characteristics as in Mihara 
(1995), while the inclusion of such feature did not im- 
prove the fit. Second, a pulse-phase spectroscopy on the 
PDS data alone shows again that a ~ 27 keV CRF is not 
detected, but only the ^ 57 keV CRF fits well the data 
(Orlandini et al. 1997b). Also a ratio between PDS spec- 
tra taken during the pulse peaks and those taken during 
the valleys did not show any variation at ~ 27 keV, but 
only at ~ 57 keV. 

But the most striking result has been obtained by per- 
forming the ratio between the PDS Vela X-1 and Crab 
count rate spectra. This ratio has the advantage of mini- 
mizing the effects due to the detector response and uncer- 
tainties in the calibrations. As we can see from the upper 
panel in Fig. 3, there is an evident change of slope at ~ 57 
keV. The same ratio performed on the Hercules X-1 spec- 
trum shows a change of slope at ~ 40 keV, just its well 
known cyclotron energy (Dal Fiume et al. 1997). 

To better enhance the effect due to CRFs, we multi- 
plied the Vela X-1 /Crab Ratio by the functional form of 
the Crab spectrum, i.e. a simple power law with index 
equal to 2.1. In this way we point out deviations of the 
Vela X-1 spectrum from its continuum without making 
any assumption on its form. The result is shown in the 
second panel of Fig. 3, where it is evident that deviations 
occur only at E ^ 35 keV. 

Finally, in the lower panel we show the ratio between 
the previous function and the Vela X-1 continuum func- 
tional, with ai,a2, and Ec taken from Table 2. In this way 
we enlarge all the effects due to line features, although we 
introduce a model dependence. We can clearly see the ab- 
sorption feature at ~ 60 keV, that we fit between 40 and 
80 keV with both a Gaussian (shown in the figure) and a 
Lorenzian. The results are shown in Table 3 and seem to 
indicate a preference for the Gaussian shape. 

Because of these considerations we will interpret the 
~ 57 keV CRF as the fundamental cyclotron resonance. 
This corresponds to a surface magnetic field strength 
of 4.9 X lO'^^ Gauss, in agreement with theoretical con- 
siderations about a high magnetic field in Vela X-1 
(Borner et al. 1987). 

At the cyclotron resonance frequency Wc, electrons at 
rest absorb photons of energy UcOc- For moving, thermal 



Table 3. Fit results on the Vela X-1 CRF. Errors represent 
90% confidence level. 





Gaussian 


Lorenzian 


^rmcyc (keV) 


60.4 ± 0.7 


59.2 ± 0.6 


CTcyc (keV) 


11.6 ±0.5 


22 ±1 


xl (dof) 


1.275 (12) 


6.662 (12) 



electrons the Doppler broadening Aud is predicted to be 
(Meszaros 1992) 

/ 2kT 
^'"'=''^[^) 

where kTg is the electron temperature, and mgC^ is the 
electron rest mass. The angle measures the direction of 
the magnetic field with respect to the line of sight. Outside 
the range lOc^Acod the cyclotron absorption coefficient de- 
cays exponentially, and other radiative processes become 
important. From Eq. 1 and the CRF parameters given in 
Table 2 we obtain a lower limit to the electron tempera- 
ture of ~ 8 keV, in fair agreement with the calculations 
of self-emitting atmospheres of Harding et al. (1984). A 
~ 27 keV CRF would corresponds to an electron temper- 
ature, assuming the same FWHM, of ~ 34 keV: a bit too 
high. On the other hand, by imposing fcTg ^ 8 keV would 
correspond to a CRF FWHM of ~ 10 keV that should be 
clearly visible in our spectra. 

Finally, the ^ 57 keV CRF agrees with the empiri- 
cal positive correlation found between the magnetic field 
strength and the spectral hardness, whether defined in 
terms of hardness ratio (Frontera & Dal Fiume 1989) or 
high energy cutoff (Makishima et al. 1990). 

5. Conclusions 

The wide-band average spectrum of Vela X-1 as observed 
by BeppoSAX appears to be quite complex, showing a pro- 
gressive change of slope: first at ^ 10 keV and then at ^ 35 
keV, just in the same interval where the Vela X-1 pulse 
profile changes from a five-peak to a two-peak shape. This 
is clearly visible in the single NFIs spectra, described by 
power laws with significantly different indices. We approx- 
imate this behaviour with a double power law modified by 
an exponential cutoff (the so called NPEX model). From 
the form of the residuals we added a cyclotron resonance 
in the fitting functional: the BeppoSAX data do not re- 
quire a CRF at ~ 27, as found in earlier works, but only 
at ~ 57 keV. A possible explanation could be an intrinsic 
variability in the cyclotron resonance, or an effect due to 
the change of slope of the continuum, not previously taken 
into account. A better description of the continuum, be- 
yond the scope of this paper is surely needed. 
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